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Abstract 

We consider the propagation of gravitational waves in six dimensions induced by sources living on 3-brancs 
in the context of a recent exact solution Q. The brane geometries are de Sitter and the bulk is a warped 
geometry supported by a positive cosmological constant as well as a 2-form flux. We show that at low energies 
ordinary gravity is reproduced, and explicitly compute the leading corrections from six dimensional effects. 
After regulating the brane we find a logarithmic dependence on the cutoff scale of brane physics even for 
modes whose frequency is much less than this energy scale. We discuss the possibility that this dependence 
can be renormalized into bulk or brane counterterms in line with effective field theory expectations. We 
discuss the inclusion of Gauss-Bonnet terms that have been used elsewhere to regulate codimension two 
branes. Wc find that such terms do not regulate codimension two branes for compact extra dimensions. 
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I. INTRODUCTION 



In this paper we discuss a simple six dimensional model, with two compact extra dimensions 
with matter living on two conical 3-branes. This model may well approximate an inflationary 
phase, or late-time dark energy dominated phase in scenarios with two large extra dimensions 
as proposed in the original scenario of Arkani-Hamed, Dimopolous and Dvali (ADD) Q, Qj. Our 
concern shall be how the matter on the brane backreacts on the bulk, and in particular the resulting 
gravitational waves that are generated. We find that at low energies, at leading order ordinary 
gravity is reproduced with the anticipated effective Newton constant. In addition there are a 
series of corrections suppressed by powers of R 2 E 2 where R c is the size of the extra dimensions. 
The magnitude of the leading order correction is sensitive to how we regulate the brane even at 
energy scales far below fi <C 1/1 where I is the width of the brane. This suggests an apparent 
model dependence in dealing with codimension two branes. We discuss the possibility that this 
dependence can be renormalized, i.e. absorbed into brane and bulk counterterms. Recently similar 
ideas have been proposed as being useful for understanding black hole physics Q|. 

The proposal of ADD resolves the hierarchy problem between the Planck and electroweak scales 
by having the fundamental Planck mass of order the electroweak scale, and using the dilution effect 
of gravity in extra dimensions to give rise to an effectively weak four-dimensional gravity. In the 
context of string theory this typically introduces a new hierarchy between the size of several large 
extra dimensions and additional small extra dimensions and/or the fundamental Planck scale, 
something which may or may not be natural depending on the details of the moduli stabilization. 
In the supersymmetric large extra dimension scenario (SLED) it has been suggested that the same 
features of the large extra dimensions combined with the calming effects of SUSY can also be used 
to ane^te the ™ logl ca> constant « pr o bta n flfl, (se e Q *. a recent ™). A 
detailed understanding of whether this is true requires amongst other things an understanding of 
how changes in the brane tensions, of the type that may arise in phase transitions on the brane, 
influence the bulk dynamics. It has proven technically challenging (| to answer this question for 
reasons that shall become apparent in the following. Here we shall focus on the simpler question 
of how the bulk dynamics is influenced by small matter perturbations on the brane. 

Although great progress has been made in recent years in understanding the dynamics of codi- 
mension one branes, higher codimension branes remain something of an enigma. This arises because 
at the level of GR, distributional sources of arbitrary codimension are typically singular , . As 
a result one needs to deal with 'thick' branes such as those described by regular defects from field 
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theory models |1C|, I 111 1 1 21 ] . An alternative approach is to modify Einstein gravity by the addition 
of higher derivative curvature terms which can potentially allow distributional sources to be well 
defined. In the case of codimension two branes one such approach that has been taken is to use 



Gauss-Bonnet terms 



0,Q,H,|lil. 



One of the main motivations for the use of these terms comes 
from string theory where they generically appear as a leading order quantum correction to gravity 
and guarantee a ghost- free action However, as will become apparent later this method of 

regularizing the branes with Gauss-Bonnet terms can only be applied for codimension two branes 
in a noncompact extra dimension, and is inconsistent for the more familiar case of compact extra 
dimensions. Some other inconsistencies have been suggested in the case of an axially symmetric 



JU 



k, where an isotropic braneworld cosmological ansatz seems to be incompatible with the model 



15 1 . 

In string theory we are interested in the dynamics of D-branes of arbitrary dimension, and 
here we find that at the level of the supergravity, the geometry describing all but the D3-branes 
are singular. In string theory we are to understand that string a' or g s corrections will 'regulate' 
the singularity. However, at first sight there seems to be a contradiction with expectations from 
low energy effective field theory (EFT), since this suggests that we need to understand the string 
scale physics that regulates the brane in order to make predictions about low energy dynamics of 
relevance to cosmology. It is precisely this aspect that we would like to explore in the following 
by asking the simpler question, how does the production of gravitational waves induced by sources 
living on the brane depend on the detailed physics that resolves the brane? 

We begin in section [H] by discussing approaches to regulating the brane geometries. Then 
in section IHII we introduce the background solution and calculate the tensor perturbations via a 
derivative expansion. This enables us to obtain order by order the modified equations for the grav- 
itational waves on the brane in section HVI We discuss the presence of the logarithmic divergences 
that arise in the thin brane limit and whether these can be canceled by counterterms on the brane 
or in the bulk. After reviewing in detail the Kaluza-Klein limit of the general solution in section 
IVl we compute in section IVTl the effects of adding a Gauss-Bonnet term in the bulk and show that 
this cannot be used to regulate the divergences near the brane. Finally in section lYlIl we conclude. 

II. DEALING WITH CODIMENSION TWO BRANES 

In this section we discuss our approach to dealing with codimension two branes. We follow 
closely the approach of ref. 0| . 
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A. Thick Branes 



One approach to dealing with codimension two objects is to regulate them by replacing them 
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191 ] . For instance one model of this would 



with a smooth stress energy source, i.e. a 'thick brane' 
be as a codimension two topological defect arising for example from a six dimensional Abelian- 
Higgs model |lo| . Figure ^ illustrates what the resulting geometry of the two extra dimensions 



would look like. The conical deficits at the poles are replaced with smooth geometries. 



FIG. 1: Six dimensional compactified spacetime with regularized singularities. 



Although the full behaviour of the solution in the polar regions may be quite complicated, 
provided the width I of the brane is much less than the scales of physical interest we anticipate that 
the bulk physics will only be sensitive to integrated effects over the brane, a perspective espoused 
in As in many works on this topic we shall restrict ourselves to axisymmetric solutions which 
is a consistent (if not entirely desirable) truncation. Given a six dimensional metric of the form 

ds 2 = N{r) 2 dr 2 + C(x, r) 2 d<p 2 + g^x, r)dx^dx u , (1) 

with branes located at r_ and r + > r_ , we define the four dimensional stress energy of the brane 
by 

(4) T M±) = 1 f r± drd^N{r)C^x)^-g^T^ (2) 

y/-g(r = r±) J r ± 

with r± = r± =F £• Here ( 6 )x^ is the regular stress energy source describing the thick brane which 
we assume vanishes outside the region \r — r±\ < e. Here e is related to the width of the brane by 
I = | JJ* N(r)dr\. The requirement that the metric is smooth as r — ► r± implies 



£(x,r±) = 0, (3) 

Nir^drCix,^) = T l, (4) 

N(r±)- x dr gitt ,(x,T±) = 0. (5) 
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The key observation of [18|] is that by integrating the Einstein equations over the brane thickness, 
one may infer an effective matching rule for the extrinsic curvature on the surface r = r± =F e 
in terms of the integrated brane stress energies. This matching rules play the role of the Israel 
junction conditions for codimension one branes. 

In particular if one looks at the equation we have 



(6) 



(6) 



Here K„ is the extrinsic curvature defined by Kg = ^jy g AC ' d r gcB- The dominant contribution is 
expected to come from the second order derivatives, and so on integrating around the brane we 
find 



2n£(r±)K»\ r=r - ± = ± K - I#«2#(±)) +0(e). 



(7) 



Similarly we have 



2ir£(r ± )K$\ 



tp\r=r± 



±7T- 



g\r=r± 
'd\r=f± 



(4) T ^(±) _ 1(4) T M(±)A +Q{e) 



(8) 



Note that the terms all vanish due to the initial assumption of the form of the metric. 
B. Regularizing as a Codimension One Brane 

The previous matching conditions amount to a statement about the extrinsic curvature defined 
on a codimension one surface in terms of the brane stress energy which is precisely the same physical 
information as if we really had a codimension one brane, localized for example at an orbifold fixed 
point as in the Randall-Sundrum model as shown in figure HI 




FIG. 2: Six dimensional compactificd spacetime with codimcnsion-one orbifold planes cutting off the singu- 
larity 
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Consequently we may regulate the codimension two brane as a codimension one brane located 
at an orbifold fixed point, where the Israel matching conditions would give 

K»\ r= r ± = ±\k (Mt^±) - ^ 8 tiQ T M{±)\ + Q(e) (9) 

K\ r=f± = ±\k ( (5) r^W - \^T^ + O(e). (10) 

Note that the factor of 1/2 arises because of the doubling of the extrinsic curvature at the fixed 
point. The five dimensional stress energy is related to the four dimensional stress energy by 

(5) T M±) = 1 ( (4) T M(±) + %( (4) r ,(±) + £ lESfA) (11) 

vr£(x,r±) V V K V-9\r±JJ 

(5) T tp(±) _ 1 ((*)J*P(.±) + (4) T r(±)\ _ / 12 \ 

v ir£(x,r±) V ^ r J 

C. Tensor perturbations 

One of the problems with the above matching rules is that they only contain partial information 
without a specification of ^TJT and ^T%. In the following we shall be interested in tensor per- 
turbations about a fixed cosmological background (gravitational waves), where tensors are defined 
with respect to four dimensional observers on the brane. In this case we have 8T£ = 5T^ = and 
so 

s (S) T j*(±) = 1 SP>T£W (13) 

ir£(x,r±) 

j(5) r ^(±) = o. (14) 

It is this simplification that will allow us to determine the bulk solution without further specification 
of the physics of the brane other than through the cutoff e. 



III. SIX DIMENSIONAL SOLUTIONS 



We use the index conventions that Greek indices are four dimensional, labeling the transverse 
directions, while small Roman indices are five dimensional, labeling the x^ and (p coordinates. 

The full six dimensional coordinates are represented by capital Roman indices: A = 0, • • • ,5. Our 

starting point is the action 

S = J d'x^l- - 2A - \F AB F AB \ , (15) 
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for gravity and a bulk form field Fab = 9aAb — 8bAa, {Ab bein g a U(l) gauge field) and 
cosmological constant A. For other six dimensional solutions see [^ol . 21] and for cosmology on 
codimension two branes, see for instance 22J, |23|. 



A. Background solution 

We consider a special case of (0), where fluxes are present in the six dimensional bulk and 
C(r, x) 2 = L 2 f(r), leading to the compactification of Q], with a metric of the form 

ds 2 = g AB dx A dx B = r l dr 2 + L 2 fd<p 2 + g^dx^dx" (16) 
g^dx^dx" = flVg^dxMx", (17) 

where q^ u is the four-dimensional metric. For the background, q^ = 7^ = a^r)^^, a(r) = 
(—Ht)^ 1 and / = 1 — -^r 2 — ^ — where /i, b are some constants and A is the six dimensional 
cosmological constant. The gauge field is of the form Aj^dx M = ^ Ld(p. H is an arbitrary 
reference scale that we have included to keep dimension. L is included so that we can normalize 
ip G [0, 2ir]. The proper size of the <p direction is then 2*K\f~fL. 

The properties of this solution are discussed more fully in [l], the main point is that provided 
A > 0, one can find solutions for which / vanishes linearly at two positive points on the real axis 
r = r±, and for arbitrary values of the parameters these points will be locally Minkowski with a 
conical deficit angle signifying the positions of two conical 3-branes. To make this clear one may 
define p± = \J ^ 2 ^^J ^- an d &± = = F^/'( r ±)) f° r which the near brane metric becomes 

ds 2 w d/4 + L 2 a 2 ± p 2 ±dtp 2 + H 2 r 2 ± q^ v dx^dx v '. (18) 

The conical deficit angle is defined through 27r(l — 5±) = 2irLa± which corresponds to two branes 
of tension 

T ± = K - 1 (l-La±), (19) 

where k is the six-dimensional Newtonian constant. For fixed bulk cosmological constant A, we are 
free to choose fi, b and L which gives us more than enough freedom to fit two branes of arbitrary 
tension. In fact we may set the flux to zero, i.e. 6 = and provided A > we can still find solutions 
with arbitrary brane tension by adjusting L and fi. The geometry of each brane is de Sitter with 
Hubble constants H± = l/r± which is only indirectly related to the brane tensions and the bulk 
flux and cosmological constant. 
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Boundary conditions 



Following the prescription of section 111 Al we now consider two codimension one branes located 
at r = r±. The normal vector to the branes is Na&x a = / _1//2 dr and the extrinsic curvature is 



V7 ro tsv f 
2V 7 / 



On the branes, the following junction condition should be satisfied: 

&[K5i-K a b ] f± =^) T « {± \ (21) 

where ( 5 )T fe a ^ is the stress-energy tensor for the gauge and matter fields introduced on the codi- 
mension one branes. The extrinsic curvature on the brane should hence satisfy 

« (5)t*(±) =+ *( (5) t »(±) _ I(5) T c(±) ™ 



Kg(f ± ) = ±- ^ ' = ±- ( Wj?=> - ^Z^fl ) . (22) 



For the background, we should therefore have: 

(5) f M±) = ± JL^TJ^ (23 ) 

(5) Tj? = ±~VW±). (24) 

B. Tensor perturbations 

We now study the cosmological perturbations around this background solution sourced by mat- 
ter on the branes. To start with, we consider only four-dimensional tensor perturbations to the 
background solution. As elsewhere we only consider axisymmetric perturbations which is a consis- 
tent truncation. We therefore have q^ u = 7^ + h fMU (x )1 ,r) with j^h^ = and 7^" h aU] ^ = 0. 

Working in this gauge, the only non-trivial contribution to the Ricci tensor is 

{%) R^ [9AB] = {4) R»u M - H 2 (3/ + rf) - 2H 2 r/V>,r " \r 2 H 2 (/V,r), r , (25) 
the perturbed Einstein equation is therefore 5^Rl = 

dr (fd r K) + - r fd r K = ^ 6R» = B h£. (26) 



The operator □ may be expressed as 



Bhft, 



□ - 2H 2 



K, (27) 



□ being the four-dimensional Laplacian: □ = j^DaDc 



Boundary conditions 

At the perturbed level, we consider some matter with stress-energy 5 (SYrff^ = g ac (f±) 5 
on each brane. Since we concentrate for now on four-dimensional tensor perturbations, we impose 
5 (5>rf ±) = 0, 5 W*) = 5 ( 5 >C (±) = 5 ( 5 >T£ (±) = and 5 = 0. 

From the junction conditions, we should have: 5K%{r±) = ±5 6 fyrfj^ with 

8K* = ^-BrK (28) 
6K% = 0, 5K* = 0, SKft = 0. (29) 



This relation implies: 



2 v ' r 



= ±^(6> r K±) ) (30) 



we recall that ^ = ^ a h av , whereas 5 ( 5 >T fe a(±) = g ac (r ± ) 5 , hence 5 ( 5 >T^ (±) = 

r/2 1 -a 7^ a <5 ^ 5 ^r a ^. From now on, we use the notation ^hjf^ = 7^ Q <5 ^^Tcw, and hence the junction 



condition is: d r hy(f±) = ± K , ( 5 V{5 . Furthermore, we may recall from the relation (1131) 

H 3 f^.y/f(r±) 

between the five and the four-dimensional stress-energy tensor that ( 5 V/f = 1 (4) T /4±) 

we hence have the boundary condition in terms of the four-dimensional stress-energy tensor: 

8 ^> - ' 4VgW - (31) 

IV. LOW-ENERGY EXPANSION 



At low-energies, we may consider that the contribution from the four-dimensional derivatives to 
be small in comparison to the r derivatives. We can hence consider an expansion in the operator 
□. The way of solving the Einstein equation (|26j). will be similar to the RS case. For that, we 
will express the solution of ()26|) as an expansion in □ and in what follows we omit any indices to 
lighten the notation and h will designate hit, and similarly for (^V^: 

Mr,z^)= W-^V 1 K(r,xn, (32) 

n>0 ^ ' 

where we consider each h n to be of order zero in \H/H 2 . So each h n have a similar weight in the 
expansion, but they are weighted by a factor (\Z\/H 2 ^ n 1 which makes their effective contribu- 
tion smaller and smaller in the low-energy regime. In particular we will consider the sum to be 
dominated by the zero mode ho(r,x^). 
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We can now solve the modified Einstein equation for each h n (r,x), n > 0: 

4 1 

d r g n (.r,x) + -g n (r,x) = --^h n -i(r,x), (33) 

with the notation /i_i(r, x) = and g n = fd r h n . 

Each mode should satisfy as well the junction conditions. Using the constraint (|3U|) . we therefore 
have the boundary conditions for each mode: 

d r h n (f ± ,x^) = Vn>0,ra/1 (34) 
a ,Mr ± ^) = ± „ LH2 : im) (35) 

At this point we should stress that in dealing with functions of □ (or □) we always have to be careful 
that there are implicitly homogeneous solutions of the equations. For instance in dealing with an 
equation of the form [□ + aD 2 ]/i = T we shall write its solution in the form h = n _1 (l — ad + 
0(D 2 ))T. Whilst this is a particular solution, we should also include the homogeneous solution 
satisfying [□ + aD 2 ]/i = 0. We shall take it as read in what follows that these homogeneous 
solutions should be included, and hence concentrate on the particular solution. We may point out 
that this method would give rise to precisely the same result as obtained in j^j] if it was applied 
to a codimension one brane. 

A. Zeroth order 

The Einstein equation for the zero mode Q33|) can be easily solved (in particular, it does not 
depend on /), and the solution is simply 

d r h Q {r,x») = ^-±. (36) 

The constant Dq may be fixed using the junction condition ()34j) which fixes Dq = 0, and so: 

ft (r,x'*) = Co(^). (37) 

As we shall see in what follows, Co will be fixed by the junction condition for the first mode since 
the zero mode acts as a source term for the first mode. 

B. First order 

The first order mode can be found by solving 

o r gi + -gi = -^ho = ^ — ' ( 38 ) 
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the solution of this equation satisfies 



f 3/r /r 4 



We work for now in the region r_ < f_ < r < f + < r+, and ^h(r,x^) is hence regular. The 
constants Cq{x^) and Di(x^) should be fixed from the junction conditions as follows. The junction 
condition for the first mode is given in l)35j) . We therefore have: 

^ = rt&jjk -fi) ( f - ^ + * <4,T ' +1 ) (40) 

f ~ (4) T H , f (4L(+) 



'KLti^ \ r_ — r 

so that 



Dl ( x n = Lff2 ; - + g3 s (f + ^ + f_ ^vw ) , (4i) 



1 (f 3 _ r 3) (4) r ,( + ) + | (f 3 _ r 3) (4) T ,(-; 



(42) 



r " vrLi/ 2 (fi - 7=3 ) f(r) 

In the low-energy limit, we are interested in the expression of the zero mode, which is finite in the 
limit where e — ► and the same on both branes. Its contribution to the gravitational waves is 
given by: 

^ V> ± ,^) = ( ( % (±) + | (4 V^ } ) (43) 

3r 2 

On the brane, the Hubble parameter is fixed to H± = l/r±, by making the choice H = H + or 
H = iT_ , we can work in terms of the proper coordinates on either of the branes. By making this 
choice, we have: 

K± = K = TTT^ TT K. (45) 

2irLH 2 (rj. - r 3 _ ) V ; 

C. Effective Newton constant 

The effective gravitational coupling constant we have derived is exactly what one expects by 
naively integrated out the action according to the usual argument 

J d^xsTg^R = J d i xdr^^-g^g^R^ + ... 

= 2ttL [ d 4 xdrr 2 H 2 ^— q^R, v + ... 
J 2k 



ttLH 2 {r% - r*_) J d 4 x^ir (4) ^ + ■ ■ ■ • 
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and so 



K± = K 



(46) 



2-kLH 2 (Hj. - r 3 ) 
In the limit where r + S> r_ , the zero mode couples uniquely 

to the matter on the brane at r = r + : 

2k 



(4) r M(+) 



□ 



(47) 
(48) 



This is similar to the single brane limit of the Randall- Sundrum scenario where conventional gravity 
is recovered on the positive tension brane only. 

D. Second order 



In order to understand the behaviour of the first mode hi, we need first to constrain the 
remaining degree of freedom in hi by imposing the boundary conditions on the second mode. 
The first mode is given by the integration of (|42|). In particular we may use the relation: 



1 



-dr 



„6— m 



r 6 f(r) 



dr 



log \r 



(49) 



t//(r<)=0 

for any integer — 1 < m < 6, where the integral has been performed by recognizing that the 
integrand is a ratio of two polynomials, with the numerator one being of lower order. Using this 
relation in the integral of (|42|). we hence have the expression for the first mode 



hi = Ci+J^Ai log 



r — n 



( r ~3 _ r 3) (4) T (+) + fj_ (f 3 _ r 3) (4) T (-) 



(50) 
(51) 



irLH 2 f'(ri) (? 3 _ - r^) 

We may now use this expression to derive the second mode, which using the equation (|33|) . is of 
the form: 



fd r h 2 = 92 = -t 



Do — I r /iidr 



1 



D 2 - -C ir 3 + z(r) 



z(r 



) = -jE' 4i (- r (f 2 + f' r + T 'i) + 



r 3 — r, 3 ) log I? 



(52) 
(53) 
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From the boundary conditions, <9 r /i2 should vanish on both branes at r±. This constraint fixes the 
constants C\ and D<i to 

' - ' + 

D * = — ^r— ^ — ' (55) 

'- ' + 

leading to the following expression for the first mode: 

^ = 3 ( " ( !- } + i: A log ]r - r,] . (56) 



We may check that in general, the first mode is not finite in the thin brane limit e — > 0. In 
particular, the dominant contribution to the first mode on the branes is logarithmically divergent 

hi(r ± ) = A±loge + O(e ) 

Putting all this together, we have the following effective four dimensional equations of motion on 
each brane 

h± _ _| (, v < ± . + !« T «) . rfg ^; (r±) , «w + 0(e »). (SB, 

E. Renormalization and EFT 

The presence of the logarithmic dependence on the cutoff suggests an inherent model dependence 
in the form of the solutions. That is, it will be necessary to specify certain features of the brane 
physics in order to determine a unique solution in the bulk. At first sight one may not be so 
concerned about this, since the dependence is only logarithmic in the cutoff. However, this is 
a feature of the fact that we are doing linearized perturbations. The logarithmic divergence of 
the metric perturbations near the brane is in fact a signal of the onset of a generic anisotropic 
Kasnerdike singularity which the conical singularity is unstable to. Consequently we only expect 
this dependence on e to get worse at higher orders. We stress that this logarithmic divergence 
is not the same as the more familiar infrared logarithmic divergence of massless scalar fields on 
two dimensional spacetimes, in fact we anticipate that for any codimension we will still find a 
logarithmic dependence for gravitational waves, where the infrared behaviour falls of as a power 
for higher codimension. 
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In EFT we are used to the idea that if we are interested in physics at a given energy scale, we 
can integrate out the modes whose masses are much larger than this energy scale, and the net effect 
is to just renormalize various counterterms in the effective action. It is tempting to apply the same 
philosophy here, the cutoff e is associated with an energy scale 1 /e at which modes which describe 
the brane itself will become excited. If we concentrate on physics at energy scales well below this, 
for instance at scales set by the size of the extra dimensions for which the higher derivative terms 
we have been discussing are still important, then it is natural to expect that all dependence on the 
scale 1/e can be renormalized or canceled by various counterterms localized either on the brane or 
in the bulk. 

Consider first the possibility that the loge may be absorbed in brane counterterms that are 
local functions of the metric, curvature invariants and the matter degrees of freedom on the brane. 
This minimal possibility seems natural given that the high energy brane physics is localized at the 
brane itself. Such terms would correspond to redefinitions of the stress energy on each brane. The 
problem is that the metrics on each brane depend on the same combination of "t^' at leading 
order and so any redefinition would give the same contribution on each brane. However in general 
the log divergence on each brane is different and so no simple renormalization of the stress energy 
on each brane will cancel the logarithmic dependence. 

The next possibility is to include in the brane action functions of the extrinsic curvature. This 
is only consistent as a boundary condition if we also include higher order derivative terms in the 
bulk that increase the order of the Cauchy problem. Adding counterterms in the bulk as well as on 
the brane can allow one to cancel the log e divergence on each brane, but closer inspection shows 
that this is typically always at the price of reintroducing it in the bulk. The key point is that as 
long as the fact that the perturbations diverge logarithmically near the brane is unaffected by the 
bulk counterterms, the logarithmic divergences will still show up at the brane. However, we have 
not as yet performed an exhaustive analysis of all the possible counterterms that could be used. In 
section IVTl we discuss an alternative approach based on using GB terms and show that also does 
not remove the loge divergences. 

Position of the brane 

This apparent difficulty at reconciling these results with the intuition from effective field theory 
may be a consequence of the fact that the boundary condition approach we have used in section Hi Al 
is not adequate. For instance, in this approach we are interpreting the metric evaluate at r = f ± 
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to be the metric of the brane. However in practice the brane is a thick object smoothed over a 
region of width e. We have made the assumption that the variation of the metric across the brane 
thickness is suppressed by e. In fact this is not necessarily the case, suppose for example we assume 
that on the region \r — r±\ < e, the metric perturbation varies as h ~ A + B{r — r±) 2 . Matching this 
form to the known form for \r — r±\ > e shows us that h(r = f±) — h(r = r±) ~ LH 2 r' 2> \f'(r±)\ • 
This variation is precisely of the same order as the first correction. The implication is that the 
coefficient of this correction will vary depending on where within the brane we choose to call the 
'brane position'. In some sense one can renormalize the loge dependence into the 'position' of the 
brane. Allowing for more complicated evolution of the metric in the brane regime will allow for 
more possibilities to absorb the log dependence. However, we are left with the same conclusion, 
we need to specify in more detail the physics of the brane in the region |r — r±\ < e in order to 
make predictions for even the leading order six dimensional corrections to the gravitational wave 
propagation. 

What seems to be lacking is a consistent way of separating the physics corresponding to different 
scales in the manner of effective field theory. It seems clear that something along these lines needs 
to be developed before a good understanding of the dynamics of codimension two branes can 
be achieved. Interesting work in this direction has been done in and extending these ideas 



to arbitrary codimension 
codimension braneworlds 



aranes seems to be crucial to capturing the essential physics of higher 
25]. 



F. Higher orders in the derivative expansion 

This derivative expansion may be continued to arbitrary high order. It is important to check 
that no further divergences are introduced so that our approximations are self consistent. Actually 
we can show that the next orders are regular in the limit e — > 0. Using the expression (|52|) with the 
relations (|54j) and (|55|) for the integration constants, we see that d r h>2 is regular everywhere. Since 
<9 r /i2 vanishes on the branes, hi must be regular at that point as well, and hence the second mode 
is regular everywhere. Furthermore, if for a given mode n > 2, h n is regular everywhere, then 
f(r)d r h n +i(r) = g n +i(r) = -x JT* r'h n (r')dr' is similarly finite everywhere. Since d r h n+ i vanishes 
on the branes, the next mode h n +i = J g/fdr is therefore finite on the brane as well and regular 
everywhere. Therefore only the first mode has a logarithmic divergence, all further modes are finite 
in the e — ► limit on the branes and this derivative expansion is hence well-defined and may be 
continued to higher orders. 
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G. Compactiflcation with one brane 



It is possible to take a one-brane limit whereby the background constants are chosen so that 
the tension of one of the branes is zero and at the perturbed level the associated stress energy 
vanishes. We find that this is a well defined limit since the metric perturbations do not diverge at 
the smooth pole. 

V. KALUZA-KLEIN LIMIT 

A. Limit of the General solution 

The general solutions we have considered so far which describe warped compactifications, con- 
tain the more familiar Kaluza-Klein compactifications as a special limit. In particular we may 
obtain solutions describing d5 4 x S 2 as follows: Redefine r = H~ l + op and specify the constants 
H, b and L as 

/, = ^AiT- 5 - -W 3 , lfe 2 if 6 = W 2 + a 2 -l, L = ^, (59) 
15 o 12 o a 

where 

-^ = 2A-9F 2 . (60) 

Then on taking the limit ff-*0we find the metric 

ds 2 = ds 2 dS , + }{ P T l dp 2 + R 2 C f(p)d<p 2 , (61) 

with f(p) = f(p)/a 2 = 1 — p 2 /R 2 . This describes a direct product of four-dimensional de Sitter 
with Hubble constant H and a two-sphere with curvature radius R c - After a suitable gauge 
transformation the gauge field is given by A = — (bH 4 R c ) pdtp corresponding to a constant flux. 
Solutions describing Minkowski 4 x S 2 and AdS 4 x S 2 are similarly obtained by taking H 2 = 
and H 2 < respectively. 

B. Kaluza-Klein solution 

In the special KK limit, the exact behavior of the different modes may be computed exactly. 
This provides a useful test on the procedure we have presented, but allows us as well to understand 
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the behavior and the nature of the divergence beyond the zeroth order. The evolution equation in 
this limit Hr — ► 1 is simply 



dp yfd p h 

d P h \±R cT e 



-ah 



.(4) r (±) 



(62) 
(63) 



7rR c f(±R c Te) 

where the four-dimensional indices //, v have been omitted to simplify the notation. The solution 
is therefore of the form: 



h = CiVmip/Rc) + C 2 Q m (p/R c ), 



(64) 



where V m and Q m are the Legendre polynomials and m = —1/2+ \/l/4 + R% □. The two integra- 
tion functions C12 can be fixed using the boundary condition Q63[) . Using these, the gravitational 
waves on each brane hr 1 are sourced by the matter perturbations in the following way: 



2vr 



' g 2R +1 '" + ll ->"- ] 1 (4V(±)4 



7T 



(4) T (T) 



sm mir 



(65) 



where H m is the m th Harmonic number: iif m = XX=i ^ 1 which for noninteger m can be defined 
by H m = Jq l i* x dx. To leading order, we therefore have: 



+- 



2^i?2Q 
nRlU 



(4) r (+) + (4) r (-) 



2^ 



2vr 



(4) r (±) + ^ _ (4) T (T) 



1 + log 



+ 



2R C 



(4) r (±) + (4) r (qF) 



(66) 



The essential features are the same as those observed in the general case. First we may emphasize 
that only the first mode diverges and its divergence is very mild since it is only logarithmic in the 
cutoff parameter e. Another important feature is that this divergence only couples to the matter 
on the specific brane and not to the matter content of the other brane. 

This represents a consistency check on the procedure used in the low-energy limit to derive each 
modes separately, since we recover the same result. We recover for both the leading and first order 
in the expansion, the results obtained in (|43|) and (|57|) corresponds precisely with the result in l)66|) 
in the limit a — > 0. 



VI. EFFECT OF GAUSS-BONNET TERMS 



Faced with the fact that higher codimension distributional sources are typically singular in 
GR, one popular approach to dealing with this is to include higher derivative terms in the bulk 
which allow distributional sources to be consistent with the equations of motion. In particular, 
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for co dimension two branes the main focus has been on Gauss-Bonnet (GB) terms. The effect 
of introducing these terms in the bulk is to modify the boundary conditions in such a way that 
the singular solutions may be discarded 



261 ], For instance a simple analysis of a codimension two 
brane in uncompactified Minkowski space shows that at the perturbed level the bulk equations of 
motion are unaffected and hence the logarithmically diverging solution is still present, however the 
boundary conditions are modified so that we may ignore this solution. In this section we point out 
that as soon as we consider codimension two branes in a compactified space this procedure fails. 
The reason is simply that if we restrict ourself to the solution which is regular at one brane, this 
solution will inevitably diverge at the other brane (or pole if no brane is present). 
Introducing a GB term in the six dimensional action: 

S G B = ^j d 6 *^ ^) R ABCD (V Rabcd _ 4 (6) R AB (,) Rab + (6 )jR 2j > (6?) 

where a is a dimensionless parameter, the presence of this term modifies both the bulk Einstein 
equation and the junction condition. In particular, the modified Einstein equation is: 

^G AB + «^ (2 K% - ^TZ GB gAB j =-Ag AB + (f a c F bc - X -F CD F CD g A B^j , (68) 

with 

11% = {6) R {6) Rab - 2 ^Rac {6) R°b ~ 2 {6) R CD {6) Racbd + ®R A DEF {&) Rbdef- (69) 

In the Kaluza-Klein limit, we consider the bulk metric of the form 

ds 2 = ds 2 ds4 + fipy'dp 2 + L 2 f(p)d<p 2 . (70) 

The limit L = R c corresponds to the previous situation where the brane tension vanishes, and the 
background geometry is smooth without any singularity. Taking L < R c corresponds to the rugby 
ball geometry corresponding to two equal tension branes at each pole. The contribution of the GB 
term in this geometry is: 

>GBj4 <nGB X A \ _ I " v 



2K^ A B -~-K™5 A B ) = \ *3 , (71) 

-12^5^ 

where the indices x, y run over the two extra dimensions: x,y = <p, p. These terms may hence be 
interpreted as a redefinition of the background parameters, and in particular, the metric (|70|) is a 
solution of the modified Einstein equation l)68|) if 

1 1 



R 2 1 + Ua^H 2 



[2A - (9 + 12a^H 2 ) H 2 } (72) 



A B dx B = - yj 1 - 3R 2 H 2 + 12a^H 2 (Rc 2 - H 2 ) p dap. (73) 
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In this limit, the form of the background solution is unaffected by the presence of the GB term, 
and in particular in the solution that describes Minkowski 4 x S 2 , ie. when H 2 = 0, the GB 
terms vanish at this order. If we concentrate on this specific solution, at the perturbed level, the 
contribution of these terms is simply: 

SU GB ^ = -^Dh^, (74) 

and all other component vanish. The equation of motion of the gravitation waves is only very 
slightly modified: 

d P (fd p h^) = - (l - 2^r) UhP v . (75) 

The gravitational waves hence behave identically as in (|64[). with only a modification of the pa- 
rameter m: m = — | + \J\ + R 2 {^ ~ 2a.yJH / R 2 ) □. Although the boundary conditions might be 
modified we argue that the addition of this term can not remove the logarithmic divergence ob- 
tained in ()65U66|) . The key point is that the form of the solution ()64j) remains completely unaffected 
by the addition of the GB term. Both Legendre polynomials have a logarithmic divergence when 
p — > ±i? c and the solution is therefore of the form: 

h(p = R c -e) - log (76) 

„ , I ' C\ sinm7r Co cos?7i7r\ e . , 

Vp = -«. + ,) - (__ + ^^J log _. (77) 

Requiring that the logarithmic divergence cancels on both branes, would fix both parameters C\ 
and Ci to zero, giving rise to the trivial solution. This is clearly an unphysical restriction, and so 
we conclude that the GB terms do not regulate the branes. 

Although we have reached this conclusion for the case of linear perturbations around the Kaluza- 
Klein solution, it should be clear that the same result will occur in general. In fact it is a funda- 
mental feature of the GB terms that they do not change the nature of the Cauchy problem (i.e. 
the differential equations remain second order) j^, as a result we are always dealing with 'two 
solutions' both of which will diverge at one or both branes even in the nonlinear case. Furthermore 
the GB terms will not significantly alter the asymptotic form of the solutions near the brane since 
they only become significant when the curvature becomes large, and for conical singularities the 
curvature remains finite and arbitrarily small up to the singularity itself. Demanding that the 
metric is conical at both branes is too restrictive a condition on the space of solutions that we will 
be left with essentially a trivial solution. This fact casts serious doubt on whether it makes sense 
to regulate codimension two branes with GB terms even in the uncompactified case, since if the 
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regulating physics is local then this case should be no different that the compactified case in the 
limit in which the extra dimensions are very large. 



VII. CONCLUSION 



In the context of an explicit six dimensional braneworld model, we have considered the effective 
equations describing the propagation of bulk gravitational waves induced by matter sources living 
on the branes. We have shown that one obtains ordinary gravity at low energies, justifying the 
boundary condition approach we have used, and have explicitly determined the leading order mod- 
ifications due to six dimensional effects. We have seen that there is an apparent model dependence 
in these corrections and have discussed the possibility of renormalizing this in brane or bulk coun- 
terterms, in accordance with expectations from effective field theory. Although it seems plausible 
that the cutoff dependence may be absorbed into counterterms localized near the brane, we find 
that it is technically difficult to do so. Our work highlights the need to develop a more consistent 
picture of how to infer bulk dynamics from the boundary conditions imposed by the brane physics 
using the ideas of effective field theory. We have shown that the use of Gauss-Bonnet terms to 
regulate codimension two branes is inconsistent when the extra dimensions are compact, suggesting 
in fact that even in the uncompactified case they represent an unphysical regularization. Since we 
focussed on tensor perturbations, much of the analysis was straightforward, but many questions 
remain such as how to deal with the boundary conditions in the scalar sector, and what happens 
to the logarithmic divergences in the nonlinear theory. Many of the features of codimension two 
branes discussed here will be present for higher codimension branes, although the codimension two 
is a specific case that can not be treated the same way as higher codimensional branes j^ . 
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